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Figure 1: Left: Hadron spectrum at a = 0.09 fm[g| together with experimental values. Inputs to fix m„^, m s 
and a in QCD are m % , iuk and mq, respectively. Right: The ratio m^/m M i as a function of m uc i at a = 0.09 
fm and L = 2.9 fm[[|. The 5 data points lowest in mass(PACS-CS) correspond to m, = 702,570,412,296 
and 156 MeV (from right to left). The vertical dashed line is the physical point. 



1. Introduction 



Nowadays lattice QCD and chiral perturbation theory rely on each other for their develop- 
ments. In lattice QCD simulations, besides statistical errors, there exist several systematic errors 
such as finite size effect and finite lattice spacing effect. The effect due to the heavier u, d quark 
masses is one of such systematic errors: Typical u, d quark masses employed in current lattice QCD 
simulations are heavier than their physical values, and one has to make extrapolations of results to 
the physical quark masses using, for example, chiral perturbation theory(ChPT). Finite size effects 
in lattice QCD simulations may also be corrected by ChPT. On the other hand, lattice QCD provide 
useful tools to check the convergence of ChPT and to determine low energy constants(LECs) of 
ChPT, by varying quark masses in simulations 

In this review, we compare results related to pion physics obtained from various lattice QCD 
simulations with predictions from ChPT. We mainly consider mass and decay constant of pion. 
In addition we briefly discuss recent results of pion form factors for full QCD simulations and 
compare both momentum and quark mass dependences of form factors with predictions by ChPT. 
In this review we will not collect and compare values of LECs from various simulations and groups, 
since an excellent review on LECs of ChPT from lattice QCD has already existed^]. 



2. Recent full QCD simulations 

Both increases of computational resources and improvements of numerical algorithms enable 
us to perform full QCD simulations with vert light quark masses. The lattice spacing a(fm), the 
lattice size L(fm), the minimum pion mass m™ n (MeV) and Lm 1 ™" for recent large scale full QCD 
simulations are listed in table |], in the upper half of which, simulations with the conventional quark 
action such as Wilson-type or staggered-type quark action are collected, while those with the chiral 
symmetric action such as the overlap or domain-wall quark action are given in the lower half. 

In Fig|T](left), one of the recent results for hadron spectra from PACS-CS collaboration^] is 
presented. This result is obtained at a = 0.09 fm and L = 2.9 fm with 2+1 flavor full QCD using 
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Table 1: (Incomplete) list of recent full QCD simulations. 



Group 


a(fm) 


L(fm) 


m™ n (MeV) 


Lmf n 


Conventional quark action 


2+1 flavors 










PACS-CS[§ 


0.09 


2.9 


160 


2.3 


NDLC0I 


>0.06 


3.3 


240 


4 


BMW| 


\ 


> 0.065 


>4.2 


190 


4 


J-Lab.[g| 


0.012 


1.5-2.9 


385 


5.7 


2 flavors 










CERN-ToV[§] 


>0.05 


1.7-1.9 


300 


2.9 


ETMC0 


>0.07 


2.1 


300 


3.2 


CLS 


i 




0.08 


2.6 


230 


3 


QCDSF[g] 


> 0.072 


2.3 


240 


2.8 


Chirally symmetric quark action 


2+1 flavors 










RBC-UKQCD[[I(J 


0.11 


2.8 


330 


4.6 


JLQCDjjll]] 


0.11 


1.8 


315 


2.8 


2 flavors 










RBC[ 


12 


] 


0.12 


2.5 


490 


6.1 


JLQCD[jT3|] 


0.12 


1.9 


290 


2.8 



the non-perturbatively 0(a) improved Wilson quark action. Hadron spectra, where masses of %, K 
and Q. are used to fix the light quark mass, the strange quark mass and the lattice spacing, agree 
reasonably well with experimental values, even though the continuum extrapolation has not been 
taken yet. The minimum pion mass in this simulation reaches 156 MeV, which is almost equal to 
the physical pion mass, though the finite size effect could be large at m n L = 2.3 of this pion mass. 
In Fig .[fright), the pion mass squared in this simulation divided by the light quark mass ( rn^/niud) 
is plotted as a function of the light quark mass. As the quark mass decreases, m 2 n /m uc i positively 
deviates from the linear behavior seen in the previous result of PACS/JLQCD collaborations at 
heavier quark masses. This deviation is expected form chiral perturbation theory(ChPT) at the 
next-to-leading order(NLO). To check a magnitude of a possible finite size effect on the spectra at 
the lightest quark mass, a new simulation at L = 5.8 fm and m % ~ 140 MeV {m % L > 4), which is 
indeed considered as a "real QCD simulation", is on-going JT4]]. 

3. Pion mass and decay constant and chiral perturbation theory 

In this section we compare quark mass dependences of pion mass and decay constant from 
recent full QCD simulations with predictions by ChPT. 

3.1 2-flavor QCD and SU(2) ChPT 

We first consider lattice QCD simulations with 2 degenerate dynamical quarks. In this case, 
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SU(2) ChPT at NLO predicts 0] 

f 2/?m„ / 2Bm,.\ 1 1 

(3.1) 



m 2 , ( 2Bm a 

— = 2B I 1 + q 



fn=f<l 



167T 2 f 2 
2Bm, 



2Bm, 



8?r 2 / 2 



ln(-^)-/,iu> 



(3.2) 



where is the pion mass, f % is the pion decay constant with the f % = 132 MeV normalization, 
B and / are low energy constants (LECs) at the leading order (LO), and £3,4 (ju) are those at NLO 
with renormalization scale \i. 

JLQCD and TWQCD collaborations have performed 2-flavor QCD simulations, using the 
overlap quark, which has an exact "lattice crural" symmetry, at a = 0.12 fm and L = 1.9 fm with 



the global topological charge being fixed to Q = 0Q13Q. The minimum pion mass in the simulation 
is m K = 290 MeV, which corresponds to m n L = 2.9. Before fitting data with ChPT formula, 1 /V 
finite size effect due to the fixing topological charge[|l^] as well as the ordinary exponential type 



finite size effect[17] have been corrected by ChPT at NLO with phenomenological values for NLO 



LECs[]18p. It is found that two finite size corrections tend to cancel each other, so that a total 
correction is small for both m K and f % . Three different expansion parameters x = 2Bm q / (87T 2 / 2 ), 
x = ml/(S7C 2 f 2 ) and % = m\ / '(87T 2 / 2 ) are employed in the ChPT fit at NLO, where m n and f % 
are measured value of pion mass and decay constant at each quark mass m q . Fig.||(Left), where 
m\jm q and f n are plotted as function of m 2 , indicates that NLO ChPT fits work reasonably well 
at the lightest 3 pion masses at m K < 450 MeV, for all 3 choices of the expansion parameter. This 
fact establishes the validity of the NLO ChPT fits at small enough pion mass (m n < 450 MeV). 
Furthermore, it is noticeable that the <^ -fit describes data beyond the fitted region. The NNLO ChPT 



fit[]T9|] with the t, is found to be reasonable for all data points at m n < 750 MeV, if a combination of 
NLO LECs appeared at NNLO is fixed to the phenomenological value. However fits show that the 
NNLO correction become 30% for m K and 70% for f % of the NLO correction at m n = 500 MeV, 
and variations of some LECs from NLO to NNLO are significant. 

The twisted mass formulation of lattice QCD is defined by a combination of the Wilson quark 
action and the twisted mass term given by m q q(x)e' 9r5 ' !:3 q(x) where q(x) = (u(x),d(x)). A remark- 
able property of this formulation is that 0(a) lattice artifacts appeared in the Wilson formulation are 
absent for physical observables at the maximal twist angle, 6r = Tt/2, where Qr is a renormalized 
twist angle. European Twisted Mass Collaborations (ETMC) have employed the twisted mass lat- 
tice QCD at the maximal twist for 2-flavor full QCD simulations, by numerically tuning Qr^h/2, 



at a = 0.087, 0.067 fm and L ~ 2.1 fm[|7|, j20] |21|]. The minimum pion mass is 310 MeV, so that 



m K L ~ 3.3. The finite volume effects have been corrected by NLO ChPT[|23|, |17fl . Fig.g(Right) 
shows na dependences of arn^/ jiileft) and a/^(right) at a = 0.087 fm, where ji. = ra^sinOj? is the 
twisted quark mass[|^]. Results indicate again that NLO ChPT fits work well at m n < 500 MeV for 
m| and f n . Fits including NNLO corrections or scaling violations have also been performed pT|], 
It is then found that the NNLO fits also works but the variation of the LEC £3 is significant from 
NLO to NNLO while a change from a = 0.087 to 0.067 fm seems to be explained by ChPT if the 
scaling violation is included at NLO. 

CERN-TorVergata groups have employed the non-perturbatively 0(a) improved Wilson quark 
action for 2-flavor full QCD simulations at a = 0.052,0.072,0.078 fm and L = 1.7 ~ 1.9 fm[0]. 
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Figure 2: Left: m^/m ? (top) and / w (bottom) from the overlap quark[|T^| as a function of m\ at a = 0.12 fm, 
together with the NLO ChPT fits with 3 different expansion parameters. Right: The same quantities from 
the twisted mass QCD, taken from Ref.[E0fl. 



They have found that the quark mass dependence of m\ at m n = 377 ~ 495 MeV is almost inde- 
pendent on the lattice spacing and is consistent with the NLO ChPT formula. 

Let me summarize the current situation of the 2-flavor QCD and SU(2) ChPT. First of all, the 
NLO ChPT describes lattice data of m 2 n and f n well at m n < 500 MeV and the expected continuum 
chiral-log is now unambiguously observed on the lattice for the first time in the unitary theory. 
The NNLO ChPT may fit data beyond this pion mass, if some NLO LECs which starts appear- 
ing at NNLO are fixed to some phenomenological values. It turns out, however, that the NNLO 
corrections seems large, in particular for f n , and that values of some NLO LECs are significantly 
affected. 

I would like to give some remarks before closing this subsection. In the 2-flavor QCD simula- 
tions mentioned above, possible finite size effects have been corrected by assuming that the ChPT 
formula for the corrections are valid. Therefore it is important to check this assumption by lat- 
tice simulations. In table||, magnitudes of finite size corrections, Rq = {O(L) — 0(°°})/0(oo) for 



O = m K or f K , from lattice QCD simulation, the NLO ChPT[g3|] and the resumed NLO ChPT[]17|], 
are compared in the case of the twisted mass QCD [22]. Although the resumed NLO ChPT formula 
is roughly consistent with lattice result, more detailed investigations are needed for the definite 
conclusion. To see whether the NNLO ChPT indeed describes lattice data well, the NNLO ChPT 
fit without using phenomenological inputs should be tested. For this purpose, simultaneous fits to 
various quantities are needed to stabilize NNLO fits. For consistency, finite size corrections should 
be included in the fitting formula, instead of correcting lattice data before the fits. Finally, from the 
theoretical point of view, inclusions of the lattice artifacts in ChPT should be made to fits data cor- 
rectly for Wilson-type quarks [24, 25] including twisted mass formulation! 26], as has been shown 
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Table 2: Finite volume corrections to m % and f n from lattice QCD simulations [Q, NLO ChPT[^3|] and 
resumed NLO ChPT formulaQ 





Rm n % 






lattice 


ChPT 


resumed 


lattice 


ChPT 


resumed 


3.0 


+6.2 


+1.8 


+4.7 


-10.7 


-7.3 


-8.9 


3.3 


+1.8 


+ 0.62 


+1.0 


-2.5 


-2.5 


-2.4 


3.5 


+1.1 


+0.8 


+1.3 


-1.8 


-3.2 


-2.9 



that it is mandatory to include such corrections in the case of the staggered-type quarks [^7|]. 

3.2 2+1 flavor QCD and ChPT 

The NLO formula of SU(3) ChPT for m n and f % are given by 

— = 2B { 1 + - + [16m;(2L 8 (At) -L 5 fa)) + 16(2m z +m,)(2L 6 ( i u) -Ufa))] } 
m i I 3 Jo" J 

fn = fo{l-2^-VK + ^[8miL 5 fa) + 8(2mi+m s )Ufa)}^ (3.3) 

Hps = l6 x2f2 yjJT )> m n = 2B o m l, m% = B (mi + m s ), = — (mi + m s ) 

where Bo and fo are LECs of SU(3) ChPT at LO, L, (/x)'s are those at NLO, m/ is the degenerate up 
and down quark mass, and m s is the strange quark mass. The validity of the NLO SU(3) ChPT for 
mass and decay constant of pion in the 2+1 flavor QCD, however, is still unclear, since the strange 
quark mass is much heavier than up and down quark masses. An alternative theoretical framework 
to describe the quark mass dependence of these quantities is the SU(2) ChPT where the strange 



quark is treated as a heavy quark. The NLO formula for SU(2) ChPT are identical to eqs.( [3.1[ ) and 
but LECs B, f, ^3,4 (/^) appeared in the formula are m s dependent in this case. 
The SU(2) ChPT fit to the 2+1 flavor lattice QCD data has been first introduced by RBC- 
UKQCD Collaborations [0], and more recent analysis can be found in Ref.[||]. RBC-UKQCD 
Collaborations have performed simulations at a = 0.11 fm with L = 2.7 fm and at a = 0.08 fm 
with L = 2.6 fm, where the minimum pion mass 330 MeV {m n L ~ 4.6) at a = 0.11 fm and 310 
MeV (m K L ~ 4.1) at a = 0.08 fm. They have employed the domain-wall quark whose additive mass 
renormalization term is small, m res a = 0.003 at a = 0. 1 1 fm and 0.007 at a = 0.08. Since the theory 
is almost chiral with this very small additive mass renormalization, the continuum ChPT at NLO 
is used for the analysis, with the replacement that my — > rhf = nif + m res . They compare SU(2) 
and SU(3) partially quenched(PQ) ChPT at NLO, in order to describe m/ dependences of mass and 
decay constant of pion, while the strange quark mass is fixed to one value in their simulations. At 
a = 0. 1 1 fm, the partially quenched data of m\ and f n are simultaneously fitted by the NLO SU(2) 
PQChPT at dynamical quark masses, (m/a,m,a) = (0.005,0.04) and (0.01,0.04). Data and fitted 
lines are shown in Fig]3|, where am 1 ^ j (m av ) and f^a are plotted as a function of m y a at ra/a = 0.005 
(m^ ea = 331 MeV). Here m™ is a mass of pseudo-scalar meson composed of two valence quarks 



6 



Pion physics on the lattice 



Sinya Aoki 




0.01 0.02 0.03 0.04 0.01 0.02 0.03 0.04 

m y m y 



Figure 3: Left: w™/m av as a function of m y at m^ ea = 331 MeV in Nf = 2+ 1 domain-wall QCD, together 
with NLO SU(2) PQChPT fit lines. Only solid symbols are included in the fit. Right: Same for f xy . Both 
figures are taken from Ref . [|l0|| . 



whose masses are m x and m y , and m av = (m x + m y )/2. The NLO SU(2) PQChPT fits data well 
at m av < 0.01(solid symbols) with x 2 /dof ~ 0.3. On the other hand, NLO SU(3) PQChPT works 
only at m av < 0.01 with #/dof ~ 0.7, which can not cover the dynamical strange quark mass, 
m s = 0.04, of this simulation. Moreover the NLO correction at m % ~ 500 MeV in the unitary 
point( mi = m x = m y ) becomes 30-40% for SU(2) and 60-70% for SU(3), which indicates that the 
NLO SU(3) ChPT is not sufficient for the strange quark. The similar conclusion that the NLO 
SU(2) ChPT works better than the NLO SU(3) ChPT holds also at a = 0.08 fm for m w a < 0.016. 

A similar conclusion that the SU(2) ChPT works much better than the SU(3) ChPT at NLO 
is also obtained for non-chiral quark action by the PACS-CS collaboration^, 29, 14], who has 
employed non-perturbatively 0(a) improved Wilson quark action at a = 0.09 fm. The absence of 
0(a) scaling violations except small contributions of 0(niag 4 ) in the simulations justifies the use 
of NLO ChPT formula of pion mass and decay constant even for this quark action[|(], ^|. The 
SU(2) ChPT fits at NLO gives % 2 /&of ~ 0.4 at m n < 500 MeV. while the SU(3) fits leads to a 
10-times larger value, # 2 /dof ~ 4. In addition they have observed that the ratio of the NLO to the 
LO corrections is much larger in SU(3) ChPT than in SU(2), suggesting worse convergence of the 
SU(3) ChPT at m n ~ 500 MeV. 

The MILC collaboration employed the rooted staggered SU(3) PQChPT formula which in- 
clude 0(a 2 ) lattice artifact, in order to fit pion mass and decay constant obtained with 2+1 flavor 
rooted staggered quarks at several lattice spacings. They have found||3|, that the NLO SU(3) fit 
fails and NNLO analytic terms (without log terms) are added to fit data at m x + m y < (0.39 ~ 0.6)m, v 
where m x ,m y are valence quark masses and m s is the dynamical strange quark mass. This also sug- 
gests the failure of the NLO SU(3) PQChPT. Indeed, in this conference [f32|], it is reported that the 
NLO rooted staggered PQChPT with the NNLO continuum PQChPT work better for the SU(2) fit 
at fixed m s ~ m^ hys and for the SU(3) fit at m s < 0.6ra]f hys where mf 1 ^ is the physical strange quark 
mass. If the NLO SU(3) ChPT works well at m s much smaller than the physical value, one can 
determine the LECs of the NLO SU(3) ChPT more accurately than now. 
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Figure 4: M^/m q (left) and f % (right) as a function of M n predicted by the NLO SU(2) ChPT formula, 
eqs. (0) and @, for JLQCD-TWQCD(blue solid), ETMC(red dotted), PACS-CS(blue dashed), RBC- 
UKQCD(red long-dashed), MILC(black dot-dashed). Green vertical lines indicate the physical point, M K = 
139 MeV. A red diamond gives f K = 132 MeV at the physical point. 



Let me summarize the current status of the ChPT for the Nf ■ = 2 + 1 lattice QCD. The NLO 
SU(3) (PQ)ChPT seems to fail at the physical strange quark mass while that of SU(2) seems to 
work at m n < 500 MeV. In the latter case, the strange quark mass dependence should be interpo- 
lated for physical predictions and LECs of the SU(2) ChPT. In order to extract LECs of the SU(3) 
ChPT, it may be better to perform Nf = 3 instead of Nf = 2 + 1 simulations, keeping masses of all 
SU(3) Nambu-Goldstone boson smaller than 500 MeV. A similar strategy is already taken for the 



staggered quark[]32j], though the use of the chirally symmetric quark action such as overlap action 
is preferable for this purpose. I would like to give one remark, made by Prof. J. Gasser during this 
conference, that the NLO SU(3) ChPT may work at m K < 500 MeV, instead of m n < 500 MeV at 
fixed m s . To check this possibility, one must tune both m/ and m s in order to keep the condition that 
mK < 500 MeV. It will be interesting to perform such simulations. 



3.3 LECs and chiral behaviors 

Instead of comparing values of LECs B, f, and i\ among various simulations, in Fig.f| we 
plot M 2 n /m q (2GeV) and f n as a function of M 2 , which is given by 



ML 



m q (2GeV) 
fn 

- 



2B(2GeV) <^ 1 + 



Mi 



16k 2 f 2 



/I 



Mi 



%% 2 f 2 



In 



M L 
m 2 



i Ml 
ln( -f 



U(m n ) 



(3.4) 



(3.5) 



where we take ju = m n = 139 MeV. Data are taken from both Nf = 2(JLQCD-TWQCD[[T|] and 
ETMC@]) and N f = 2 + 1(MILC[|], RBC-UKQD[|l|] and PACS-CS|@1) simulations. Although 
results among different groups more or less converge for f n , sizable differences are observed for 
M 2 /m q . One of reasons for these differences seems to come from the overall renormalization factor 
of m q : MILC and PACS-CS collaborations, who employ perturbative renormalization factors, tend 
to give larger values of the ratio than those from JLQCD-TWQCD, ETMC and RBC-UKQCD 
collaborations, who use non-perturbative estimates for the renormalization. 
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Table 3: (Recent full QCD calculations for pion form factors. Here Q 2 = —q 2 . 



Group 


N f 


quarks 


a(fm) 


L(fm) 


£ 2 (GeV 2 ) 


«%(MeV) 


Fx 


QCDSF-UKQCD[|4| 
RBC-UKQCD[35] 
ETMC@ 
JLQCD-TWQCD[|7j 


2 
2+1 
2 
2 


clover 
DW 
TM 

overlap 


0.07-0.12 

0.11 
0.07-0.09 

0.12 


1.4-2.0 

2.8 
2.2-2.9 

1.9 


0.31 -4.3 
0.013-0.258 
0.05-0.8 
0.252-1.7 


400-1000 

330 
260-580 
290-750 


V 
V 
V 

v,s 



4. Pion form factors 

In this section, we consider form factors of the pion, where chiral perturbation theory plays 
an important role to make extrapolations in terms of the momentum transfer as well as the quark 
mass. The vector and scalar form factors are defined by 

W)\V ll \%{p)) = {p + p') tl F v {q 2 ), (K(p')\S\K(p))=F s (q 2 ), q 2 = (p - p'f (4.1) 



<>*)* 



dF x ( 



-2\ 



dq 2 



d 2 F x (q 2 ) 



, X = V,S (4.2) 

q 2 =0 



where (S) are the vector current (scalar density) in QCD, and (r 2 )v ( {r 2 )s ) is called the 
charge(scalar) radius. Recent full QCD calculations for the pion form factors are summarized 
in table |5[ 

It has been found by all groups that q 2 dependence of Fy(q 2 ) at fixed quark mass is well 
described by the pole ansatz: 

1 



in particular at small q 2 , as seen in Fig.||(Left). Moreover the value of M po i e obtained by the fit is 
closed to the vector meson mass M p at this quark mass (Vector Meson Dominance). The single 
pole ansatz leads to 

(r^J-, cy ^{r 2 ) v ^6^. (4.4) 

pole pole V / 

On the other hand, the SU(2) ChPT at NLO gives 

= +«W + bK/M , l}, <r = 55^"^. (4-5) 

which, together with the relation from the single pole ansatz, predicts 

(r 2 ) v ~A —L- ~ .22fm 2 , (4.6) 
V 5 A%fm % 

which is far below the experimental value, (r 2 )y Xp PDG = 0.452(11) fm 2 . This discrepancy mainly 
comes from the fact that the NLO ChPT does not reproduce behavior of lattice data and a single 
pole ansatz. 
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Figure 5: (Left) The vector form factor of pion Fy(q 2 ) as a function of Q 2 = —q 2 , together with the single 
pole fit (solid line) in Nf = 2 + 1 flavor domain-wall QCD at a = 0.11 fm and m 2 n = 330 MeV[||]. The 
dotted line is the single pole fit of Ref.[|^| at the same pion mass but in Nf — 2 QCD with clover quarks. 
(Right) The NLO SU(2) ChPT fit at m 2 = 330 MeV (dash-dotted line), 139.57 MeV (Solid line) and with 
the charged radius being fixed to the PDG world average(dashed line) for small Q 2 , together with lattice data 



at m n — 330 MeV(triangles) and experimental data (diamonds). Both figures are taken from Ref.[35 1 



Because of this problem, QCDSF-UKQCD collaborations! 34] give up the ChPT fit and em- 
ploy the single pole ansatz for the q 2 dependence of Fv(q 2 ). After the chiral extrapolation to the 
physical point by the form that M 2 ole = c + cim 2 , they obtain (r 2 ) v = 0.441 (19) (56) (-29) fm 2 , 
which agrees well with the experimental value. The form factor Fy(q 2 ) itself, the single pole ansatz 
extrapolated to the physical pion mass, reproduces experimental data well. 

Regardless of the success of the single pole fit, the situation is not satisfactory from the theo- 
retical point of view, since the NLO ChPT is incompatible with it. 



RBC-UKQCD collaborations] 35] have fitted Fy(q 2 ) at very small q 2 , instead of {r 2 )y and cy, 



by the NLO SU(2) ChPT formula given as 



F v (q 2 ) = l + -^[-2f 6 ^)q 2 +4^(ml,q 2 ^ 2 ) 



J4f(x,y,z) 



xH{y/x) 

327T 2 1927T 2 



y x 

In-, H x) 



(4.7) 



4 5* x-4 / ^(x-4)/x+\ 



Using fa = 0.0665(47) from the fit of the decay constant as an input, one unknown LEC £ r 6 can 
be extracted form one point at q 2 = —0.013 GeV 2 as i\{m p ) = —0.093(10), though no degree of 
freedom is left. This value of the LEC leads to {r 2 )y = 0.354(31) fm 2 at the simulation point where 
m % = 330 MeV, which becomes (r 2 ) v = 0.418(38) fm 2 at the physical point, m % = 139 MeV. The 
form factors Fy(q 2 ) reconstructed by the NLO SU(2) ChPT formula with the i\ are compared with 
lattice data at m n = 330 MeV and with experimental values at m n = 130 MeV in Fig. |5[ which 
shows reasonable agreements in both cases. 

The ETM collaboration^] has performed a little more advanced analysis, employing the 
NNLO SU(2) ChPT formula for the form factor. Using the experimental value of (r 2 )s = 0.61(4) 
fm 2 to fix some LECs through the NNLO formula, they have simultaneously fitted m n , f K and 
F v (q 2 ) by the NNLO SU(2) ChPT formula at 0.05 < Q 2 < 0.8 ( GeV 2 ), which predicts (r 2 ) v = 
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Figure 6: Left: The charge radius (r 2 )y(solid squires) as a function of M 2 , together with the NNLO SU(2) 
ChPT fit(solid line), its NLO (dot-dashed line) and NNLO (dashed line) contributions. The open circle is 
the experimental value extracted by ChPT and the asterisk is the PDG value. Right: The same for cy. 

0.438(29) fm 2 at m % = 139 MeV. An agreement in F v (q 2 ) between ChPT fit results at NNLO and 
experimental data is reasonably good. 

JLQCD-TWQCD collaborations^] have used a hybrid method for the fit: q 2 dependences 
of Fy and F$ are fitted by the single pole plus polynomials in q 2 , from which {r 2 )y, cy and (r 2 )s 
are extracted. On the other hand, as shown in Fig.^, m\ dependences of these three quantities 
are fitted by the NNLO SU(2) ChPT formula, which give (r 2 ) v = 0.409(23) fm 2 , c v = 3.22(17) 
GeV~ 4 (~ 0.0049 fm 4 ), and (r 2 ) s = 0.617(79) fm 2 . It is interesting to note that cy turns out to be 
very close to a value predicted from this value of {r 2 )y by the single pole ansatz, cy ~ 0.0047 fm 4 , 
though the convergence of the expansion is questionable from the fact that Cy LO < Cy NLO ( Note 
that Cy° = 0. ). It remains an important challenge to fit both q 2 and rr? % dependence of Fx{q 2 ) by 
the NNLO SU(2) ChPT formula. 

5. Summary 

Let me conclude this review. Now chiral-logs of the NLO ChPT for mass and decay constant 
of pion are clearly seen in lattice QCD even for data in unitary theories. Recent analysis suggest 
that the NLO SU(2) ChPT works for these quantities at pion mass less than 500 MeV, while an 
applicability of the NLO SU(3) ChPT to the physical strange quark mass seems questionable, so 
that the NLO SU(2) ChPT formula is often used even for data in 2+1 flavor QCD simulations. As 
far as pion form factors are concerned, more detailed investigations are needed for a comparison 
between numerical data and ChPT, in particular, for the convergence of the chiral expansion. In 
future more complicated quantities from lattice QCD, such as the %% scattering phase shift, should 
be compared with the ChPT predictions. 
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